Entanglement evolution in finite dimensions 
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We provide a relation which describes how the entanglement of two d-level systems evolves as 
either system undergoes an arbitrary physical process. The dynamics of the entanglement turns out 
to be of a simple form, and is fully captured by a single quantity. 



PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Yz. 

When precisely studying and manipulating the quan- 
tum world, the objects of interest such as photons [T], 
atoms 0, ions [3] or quantum dots [I] are usually only 
few in number. But quantum mechanics does allow also 
meso- or macroscopic systems to exhibit genuine quan- 
tum features such as interference [5J and classically un- 
achievable correlations [6J - the latter often lumped to- 
gether under the label "entanglement" . The description 
of these larger quantum systems, such as Bose-Einstein 
condensates, which consist of thousands of quantum ob- 
jects, requires to resort to effective properties since the 
quantum system's exact state surpasses what can analyt- 
ically or numerically be coped with. Therefore, we need 
efficient theoretical and experimental tools to character- 
ize and probe these quantum properties in terms of few 
and robust quantities [71 [HI E]- 

Furthermore, since the transition from microscopic to 
macroscopic scales entails a rapidly increasing density of 
states, and, thus, a strongly enhanced fragility of generic 
quantum features to perturbations, we have to account 
for decoherence, or rather for the time scales on which 
quantum properties can prevail, in the presence thereof. 
Thus, we seek an efficient dynamical description of the 
entanglement of open quantum systems [10 | H H H2], 

However, while sytematic progress has been achieved 
in describing the interference of ever larger quantum ob- 
jects [5], in ever more complex environments, our sys- 
tematic understanding of entanglement dynamics is still 
in its infancy. This lack of understanding is rooted in 
entanglement being a nonlinear function of the system's 
density matrix. Thus, coherence is a necessary though 
by no means sufficient criterium for entanglement to pre- 
vail. Consequently, the only systematic results on entan- 
glement dynamics were hitherto available for the smallest 
- microscopic - quantum systems which can harbour en- 
tanglement, i.e. for pairs of qubits |13j . Here, we present 
a first general result which describes the entanglement 
dynamics of two-party quantum systems with arbitrary 
(finite) dimensions of their components. 

The setup we consider consists of two initially entan- 
gled d-level systems, of which one undergoes an arbitrary 
physical process, in general some open system dynamics 
in which it interacts with uncontrolled and not measur- 
able degrees of freedom of its environment. Such pro- 



cesses are often referred to as channels, maps or super- 
operators [HI [T5J [IB]. We denote them by $. Starting 
with a pure state |x) the system's final state then takes 
the form 

p' = (l®$)| X )( X |, (1) 

which now needs to be characterized in terms of its entan- 
glement content. As compared to the simple case d = 2 
in |13j . the following complication arrises: p' may be a 
mixture of pure states \ipi) which are living on differ- 
ent, strict subspaces of the d-dimensional space, i.e. the 
Schmidt rank of p' may drop below d, but this no longer 
implies separability. We will here focus on the entangle- 
ment evolution on time scales while p' preserves its initial 
Schmidt rank, though will see that the infered dynamics 
has some bearing also for the entanglement evolution on 
longer time scales. 

We quantify the entanglement exhibited by p 1 using 
G-concurrence [17] . which reduces to concurrence [TS] 
when restricting to two two- level systems. For a pure 
state |x) it is the geometric mean of its d Schmidt co- 
efficients Xi {i — 1, ...,d). However, when \\) is given 
as \X) = J2i j=i Aij\i)\j) with basis states \i) and |j) for 
the respective subsystems, G-concurrence is more conve- 
niently evaluated by 

G d (\ X )) = d[det(A^A)] 1/d . (2) 

For mixed states p, G-concurrence is calculated through 
the usual minimization procedure of the ensemble's aver- 
age, Gd(p) = inf Y^iPiGd{\4>i))i over all possible decom- 
positions into pure states, i.e. p — J2iPi\ ( f > i}(4 > i\ (with 
Pi > and ^2iPi = 1). An attempt for an analyti- 
cal result of this optimization procedure similar to the 
one done by Wootters for concurrence |18j yielded com- 
putable upper and lower bounds |19j . The same frame- 
work also allows to consider (d x /)-systems, d < /, since 
pure states (and the pure state vectors that form a de- 
composition of a mixed state) of such a system can ex- 
hibit entanglement in at most d levels, by virtue of their 
Schmidt- representation . 

Exploiting some specific algebraic properties of G- 
concurrence, we can largely follow the line of argu- 
ment when describing entanglement dynamics of 2 x 2 
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states [T5] . In order to evaluate the entanglement of the 
final state p' , we first need to express the initial pure 
state |x) as the result of a so-called filtering operation 
M x [5D] acting on either party of a maximally entangled 
state |<&), i.e. 

\X) = (M X ®1)\*). (3) 

M x = Vd>Y2ij=i A-ij here acts on the first subsys- 
tem of |<f>) = Yln=i \n)\n)/Vd. Note that such filtering 
is always possible, for arbitrary |x), given its matrix rep- 
resentation Aij . The channel represented by the filtering 
operation M x and the state |x) are isomorphic [2Tj . 

Given the filtering M x and the channel $, which act 
on different and possibly spatially separated parts of the 
system, the temporal order of their execution must not 
be of influence to the final state, and we can exchange 
their order in our representation of p' , 

p' = (M x (g>t)p $ (Ml®l), (4) 

where we introduced the result of the channel $ acting 
on the maximally entangled state, p$ = (1 ® $)|<&)(<i>|, 
which is mixed in general. Again, state p$ and chan- 
nel $ are related to each other via the Jamiolkowski iso- 
morphism [21J . Summarizing these two steps above, we 
have transformed the initial setup via the Jamiolkowski 
isomorphism into the dual one, where the role of states 
(p$ replaces \x)) and channels (M x replaces $) is inter- 
changed. As in [T3] one may also arrive at this point by 
inserting an intermediate teleportation procedure before 
the d-level system undergoes the action of the channel. 

Having arrived at this particular form Q of the 
final state p', our chosen entanglement measure, G- 
concurrence for a (d x d)-system, factorizes much as 
concurrence does for the (2 x 2)-case [T3]. That is, 
G-concurrence exhibits the particular property that a 
single operator acting on either one of the subsys- 
tems simply factors out \17\ [22]: G d [(M ® l)\ip)] = 
\det(M)\ 2/d G d (\^)). Applying this to the dual form of 
the final state Q (in this very form [23]), and realizing 
that the determinant of M x relates to the G-concurrence 
of the initial state \x), yields our core result 

G4(l®$)| X )(x|]=G d (|x))G4(l®$)|$)($|] . (5) 

The entanglement of two d-level systems in terms of G- 
concurrence evolves equally for all pure states |x), is 
solely given by the evolution of a maximally entangled 
state, and merely rescaled by the initial entanglement. 
This effectively reduces the vast space of initial condi- 
tions to a single one. 

For general, i.e. mixed, initial states po and/or two 
one-sided channels $i (g> $2, the convexity property of en- 
tanglement monotones [T7l |2"41 [23] leads to an inequality 



instead, which provides an upper bound 

G rf [($i®$ 2 )po] <G d ( Po )x 

x G d [($!<8> l)|$}($|]G d [(l®$ a )|$}<$|] . (6) 

Here equality holds, for example, in case of pure initial 
states and two channels, if the filtering operation M x and 
either of the channels, $1 or $2, commute - such that the 
order of the execution of the channels and of the filtering 
can be interchanged, in order to achieve a factorization 
as found above for one-sided channels. 

Let us contemplate the physical content of eq. (|5j: 
Since G-concurrence - computed from the product of all 
Schmidt-coefficients of a state - vanishes if at least one 
of the Schmidt coefficients is zero, it measures entangle- 
ment in exactly d levels. Whenever the bipartite system 
under study undergoes some dynamics which induce a 
redistribution of amplitudes to strictly less than d levels, 
G-concurrence vanishes, while entanglement on a strict 
subset of levels may prevail. For a channel with a con- 
tinuous time evolution, e.g. governed by a Lindblad-type 
equation, the system's state evolves continuously within 
the set of states. Thus, when initially prepared in a state 
with non- vanishing G-concurrence, it is guaranteed to ex- 
hibit non- vanishing G-concurrence at least during an ini- 
tial, finite time interval, and the time evolution of entan- 
glement during this period is described by our result (|5| . 
At later times, one has to rely on a hierarchy of entan- 
glement monotones Ck(\ip))i k < d, for example those 
which were defined along with G-concurrence (being the 
last member of the hierarchy G d = C d ) [17]. Similarly 
to G-concurrence, they are computed from the sum of all 
different products of k Schmidt coefficients, and hence 
capture the entanglement of exactly k levels. For them 
we similarly derive 

C k [(l®$)\ X ){x\] <C fc (| X »C7 fc [(l®$)|*)<$|] (f) Vfe , 

(7) 

with an additional binomial coefficient |26j which ac- 
counts for the different possibilities to select k out of 
d levels. 

Finally, we conclude with a remarkable observation 
which suggests that the time evolution of G-concurrence 
is intimately related to the evolution of the entanglement 
as quantified by concurrence [37]: For this we choose an 
entangled pure state with non-vanishing G-concurrence 
and expose one of the cHevel systems for a time t to 
a depolarizing environment. Such environments com- 
pletely destroy the information about a system's state 
without any bias and thus turn it into a totally mixed 
state. When viewing the quantum system as an infor- 
mation carrier, this means that all possible independent 
errors occur with the same probability. The dynam- 
ics of G-concurrence during this procedure is then, by 
virtue of (|5| , entirely determined by the dynamics of the 
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Isotropic states (d=5) 




FIG. 1: Evolution of concurrence (solid) and Schmidt num- 
ber (dashed) as one party of the maximally entangled state 
|$) = J2n=i \ n )\ n ) I 'Vd, (d = 5) passes through a depolariz- 
ing channel with rate V (resulting in an isotropic state). The 
initial decay of entanglement is well approximated (linearly) 
by the entanglement at time (C(0)) and when the Schmidt 
number drops for the first time (here from 5 to 4, with the 
initial entanglement reduced by a fraction l/(d— 1), from 1 
to 3/4), which is exactly when G-concurrence vanishes. 

maximally entangled state. Depolarizing one subsystem 
thereof produces isotropic states [25], which can be pa- 
rameterized by 

^ = ^^(1_|$)<$|)+F|$)($|. (8) 

The state evolves along the line segment parametrized 
by F according to the solution of the Lindblad equation 
(Markovian environment assumed) with interaction rate 
r, what results in F(t) = [1 + (d 2 - 1) exp(-2dTt)]/d 2 . 
Thus, initially with F(0) = 1 the maximally entangled 
state is recovered, and thereafter the state closes in to- 
wards the totally mixed state l/d 2 , at F = l/d 2 , asymp- 
totically in time. The condition that an isotropic state be 
of Schmidt number k [3pJ [30] , namely if and only if k— 1 < 
Fd < k [30J , then determines the points in time when the 
Schmidt number drops. Thus, the drop from Schmidt 
number k to k — 1 occurs at time tf. = ln[(d 2 — l)/(dk — 
d— l)]/(2c?r) and hence entanglement measures quanti- 
fying only entanglement in k levels vanish. In particular, 
for F(t 2 ) = l/d with t 2 = ln(d + l)/(2dT), the state 
turns separable, whereas for F(td) = (ef— l)/d, at time 
t d = ln[(d 2 - l)/(d 2 -d- l)]/(2dT), G-concurrence dis- 
appears. The inverse of this depletion time of G- 
concurrence then constitutes a characteristic decay rate 
Tq compared to the decay rate Tc of the generalized con- 
currence C [27] (which vanishes if and only if the state 
is separable). Due to the high symmetry of isotropic 
states the infimum optimization for the concurrence of 
these high-dimensional mixed states can be carried out 
analytically |31j , and, when normalized to the initial con- 



currence, yields 

C( PF ) = (Fd-l)/(d-l), (9) 

for l/d < F < 1, and C(pp) = otherwise. For the 
chosen dynamics, that are encoded in F(t), it decays 
exponentially with a small offset, that guarantees it to 
vanish at time t 2 . Its time derivative then determines 
the rate Tc(t) at which concurrence decays. Comparing 
the characteristic decay rate Tq of G-concurrence to the 
decay rate Tc(t) of concurrence, at short times as illus- 
trated in Fig. [T] yields a proportionality factor which for 
any initial pure state only depends on the system dimen- 
sion: 

r G _ d d>i 

r c( o) = ((i+1)lo (^ ~ • 

This indicates that the dynamics of G-concurrence and 
concurrence are much interrelated with one-another al- 
though G-concurrence detects only a very specific type of 
entanglement, that is entanglement in exactly d levels. It 
also suggests that G-concurrence already exhibits most 
properties of the evolution of entanglement in general - 
in particular at short times, when entanglement control 
is most crucial for possible applications. 
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